Magnetic behavior of curium dioxide with non-magnetic ground state 
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In order to understand magnetic behavior observed in Cm02 with non-magnetic ground state, we 
numerically evaluate magnetic susceptibility on the basis of a seven-orbital Anderson model with 
spin-orbit coupling. Naively we do not expect magnetic behavior in Cm02, since Cm is considered 
to be tetravalent ion with six 5/ electrons and the ground state is characterized by J=0, where J 
is total angular momentum. However, there exists magnetic excited state and the excitation energy 
is smaller than the value of the Lande interval rule due to the effect of crystalline electric field 
potential. Then, we open a way to explain magnetic behavior in Cm02. 

PACS numbers: 75.20.Hr, 75.40.Cx, 71.70.Ch, 71.27.+a 



Curium is one of transuranium elements and it takes 
trivalent or tetravalent ion state. It is easy to understand 
that local / electron number n is equal to 6 and 7 for 
Cm^"*" and Cm''+ ions, respectively. By following the 
well-known Hund's rules, we obtain that the ground-state 
multiplet for Cm"^"*" is characterized by S—7/2 and L=0, 
while for Cm''"'", the multiplet is expressed by S'=3 and 
L=3, where S and L denote spin and angular momenta, 
respectively. When we further include the effect of a 
spin-orbit interaction, the ground state for Cm'^"*" with 
n=7 is characterized by J—S—ljl^ which is the same as 
that for Gd'^"'" ion. On the other hand, for the case of 
Cm'*"'" ion with n=6, the ground state is just the singlet 
with J=0. Thus, the magnetic behavior should be quite 
different between Cm'^+ and Cm'*"'" ions. 

Now we pay our attention to experimental results on 
Cm02)^ which are in a puzzling situation. When we 
naively deduce the valence of Cm ion from our experi- 
ence in other actinide dioxides. Cm should be tetravalent, 
leading to non-magnetic ground state. However, from 
the neutron diffraction experiment, the effective moment 
Mcff ~ 3.4/LtB has been observed, where /xb is the Bohr 
magneton. In accordance with the large value of /ioff, 
it has been also observed that the magnetic susceptibil- 
ity follows the so-called Curie- Weiss behavior. At the 
first glance, these experimental results seem to contra- 
dict non-magnetic ground state of Cm'*"'" ion. 

A possibility to resolve the contradiction is to consider 
the effect of Cm'^"^ magnetic impurity. However, the 
amount of Cm'^"'" deduced from the quantity of oxygen 
is apparently smaller than the value to explain the mag- 
netic moment. Furthermore, there is no evidence in the 
neutron diffraction for the superlattice peak correspond- 
ing to a long-range rearrangement of oxygen sublattice to 
accommodate a mixture of Cm"^"'" and Cm*+ ions. From 
a theoretical viewpoint, it is an interesting issue to con- 
sider a mechanism to obtain magnetic behavior from the 
non-magnetic ground state. 

One way to understand the problem is to consider seri- 
ously the magnetic excited state of Cm'*"'" ion, as already 
mentioned in the experimental paper li Namely, when the 
magnetic excited state exists just above the non-magnetic 



ground state with a small excitation energy, it may be 
possible to understand magnetic behavior even for the 
non-magnetic ground state. This scenario seems to be 
quite simple, but we believe that it is worth to examine 
faithfully the scenario from a theoretical viewpoint. 

In this paper, first we analyze the local /-electron state 
of Cm*+ ion on the basis of the Hamiltonian including 
Coulomb interactions, spin-orbit coupling A, and crys- 
talline electric field (CEF) potential. If we ignore the 
CEF potential, the excitation energy IS.E is given by 
^E=\l% in an LS coupling region, while it is expressed 
by lS.E=7\/2+rj in a j-j coupling region, where rj is the 
energy shift depending on the CEF potential. When we 
further include the CEF potential, we find the relation of 
AiScxA^ in the LS coupling region, which is smaller than 
the value in the Lande interval rule. Then, we evaluate 
the magnetic susceptibility on the basis of a seven-orbital 
Anderson model by exploiting a numerical renormaliza- 
tion group technique. It is found that the magnetic be- 
havior actually appears for a realistic value of A even if 
we assume Cm''"'" state. 

Let us first discuss the local /-electron state. The 
model is given by^ 



Hi — ^ ^(-Bm,m"^(T,cr' + A^; 



m.(7,m' ,(T' 



) fma fm' c 



m.rn cr.a 



mi ~7n4 a, a' 

where Bm.m' is the CEF potential, a=-\-\ (—1) for up 
(down) spin, fma is the annihilation operator for / elec- 
tron with spin a and ^-component m of angular momen- 
tum £=3, 5ccr' is the Kronecker's delta, A is the spin- 

Or bit coupling, C m.±l,m.=bl = ±m/2, Cm±l,Tl,m,±l = 

■\/12 — m(m ± l)/2, and zero for the other cases. The 
Coulomb integral / is expressed by the combination of 
Slater-Condon parameters such as F^, F'*, and 

Here we briefly explain the parameters of the local 
Hamiltonian Hf. Concerning Slater-Condon parameters, 
first we set F^^^IO eV by hand, since we are not in- 
terested in the determination of the absolute value of 
the ground state energy. Others are determined so as 
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FIG. 1: Excitation energy AE vs. spin-orbit coupling A for 
n=6. Thick solid curves denote the calculated results with 
different values of the CEF potential W. Broken lines indi- 
cate the two limiting situations as AE=X/6 for small A and 
AE—7\/2 + ri for large A, where r; is the energy shift depend- 
ing on the CEF potential. 



to reproduce excitation spectra of U^^ ion with two 5/ 
electronsj^ Here we show only the results: F^=6.36 eV, 
i^4=5 53 and F^=4.13 eV. As for spin-orbit coupling 
A, it may be possible to use an experimental value, but 
in this paper, we change the value of A as a parameter 
to control the situation from the LS coupling to j-j cou- 
pling schemes. Since the fluorite structure belongs to Oh 
point group, Bm,m' is given by using a couple of CEF pa- 
rameters, and Bg, which are traditionally expressed 
as Bl=Wx/l^ and Bg=VF(l - |a;|)/180, respectively.^-^ 
Here W determines the energy scale of the CEF poten- 
tial, while X specifies the CEF scheme. 

Now we summarize the experimental facts on the CEF 
energy levels of actinide dioxides: For UO2, the ground 
state is triplet and the first excited state is F;^ dou- 
blet with the excitation energy 150 meVJ For Np02, the 
ground and first excited states are, respectively, Fg 
and Fg quartets with the excitation energy 55 meV;^ 
For Pu02, the ground state is F^ singlet, while the first 
excited state is F4' triplet with the excitation energy 
123 meVj^ii^ In order to reproduce totally the above 
CEF level schemes, we have chosen PF=— 10.5 meV and 
x=0.62 in the discussion of the ground state of Am02.^^ 

In Fig. 1, we show the excitation energy AE vs. A 
for various values of W . Note that at A=0, the ground 
state is highly degenerate. Here we set AE=Q at A=0 in 
order to visualize how the degeneracy is lifted by the in- 
clusion of A. First we consider the case of W—Q, in which 
the ground state is characterized by J=0, while the first 
excited state is given by J—1. In this case, we obtain 
AE=X/6 in the LS coupling region from the so-called 
Lande interval rule. On the other hand, in the j-j cou- 
pling region, the ground state is expressed by the closed 
shell of j=5/2 sextet and the exited state is given by 



the one-electron excitation from j=5/2 sextet to j=7/2 
octet, leading to the J=l excited state. Thus, we obtain 
AE=7X/2+ri in the j-j coupling scheme, where rj denotes 
the energy shift depending on the CEF potentials. 

Now we include the effect of the CEF potential. As 
shown in Fig. 1, the excitation energy becomes small 
when we set W=—10.5 meV, which has been determined 
so as to reproduce the CEF energy levels of other actinide 
dioxides. Both the ground state and the first excited 
state energies are decreased by the CEF potential, but 
the reduction of the first excited state energy is large in 
comparison with that of the ground state energy. From 
the analysis of our numerical results, we find the relation 
of AE oc A"*, suggesting that AE is suppressed for small 
A in comparison with the relation of AE=X/6 for 1^=0. 
When we further increase the magnitude of W, the exci- 
tation energy is getting smaller and smaller, as observed 
in Fig. 1. Thus, it is concluded that the excitation energy 
of the electron system is smaller than we have naively 
expected due to the combined effect of Coulomb interac- 
tions, spin-orbit interaction, and the CEF potential. 

Let us now include the hybridization between localized 
and conduction electrons. The Hamiltonian is the seven- 
orbital Anderson modelfi^ given by 

^ = I]ekcLck. + + h.c.) + Hu (2) 

k,a" k,a",m 

where et denotes conduction electron dispersion, Cko- in- 
dicates the annihilation operator for conduction electron 
with momentum k and spin cr, and Vm is the hybridiza- 
tion between conduction and / electrons. 

Note that we consider only Ou single conduction band 
with xyz symmetry composed of oxygen 2p electrons. 
Since oxygen ions surrounding actinide ions are located 
in the [1, 1, 1] direction, there should exist a conduc- 
tion band composed of 2p electrons with xyz symmetry. 
This picture seems to be consistent with band-structure 
calculation^ but we assume the ignorance of and 
t2n bands. Here we note that the hybridization occurs 
between the states with the same symmetry of local /- 
electron state. Since the conduction band has xyz 
symmetry, we set 'Vi=—V-i=V and zero for other m, 
where V is fixed as V^=0.05 eV. In order to adjust the 
local /-electron number n, we appropriately change the 
chemical potential in the actual calculation, although we 
do not explicitly show such a term. A half of the band- 
width of Cu conduction band is set as 1 eV, which is the 
energy unit in the following calculations. 

In this paper, we analyze the seven-orbital Ander- 
son model by using the numerical renormalization group 
(NRG) technique J^ii^ In the NRG calculations, in order 
to consider efficiently the conduction electrons near the 
Fermi energy, the momentum space is logarithmically dis- 
cretized and the conduction electron states are character- 
ized by "shell" labeled by iV. The shell of A^=0 denotes 
an impurity site including / electrons. The Hamiltonian 



3 



10^ 



3 
o 

C 1 

mo 

c 
o 

w 
O 



10' 



: W=-105meV 


1 1 


A.=0 J 1 

* — ►^1=0.2 / / 

1^1=0.3 f / 
.* — *X=OA 1 1 1 

1 1 1 1 ™ 1 1 ' 1 Ll 


/ 1 



0.020 



10^° 10^ ld« ld' ld' Id' Id' ld' ld' ld^ 10° 

T(eV) 



FIG. 2: Curie constant C vs. temperature T for various values 
of A. 
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FIG. 3: Inverse susceptibility 1/x vs. temperature T for var- 
ious values of A. The CEF potential is set as W=— 105 meV. 



is transformed into the recursion form as 

Hn+I = VA-Hn + tN y^(cjv^CAf+lCT + cjv+lcrCjVcr), (3) 
a 

where A is a parameter for logarithmic discretization, 
CNa denotes the annihilation operator of conduction elec- 
tron in the iV-shell, and tj^ is "hopping" of electron be- 
tween N- and {N + l)-shells, given by 



tN = 



(l + A-i)(l-A-^-i) 



(4) 



2^/(1- A-2A^-l)(l-A-2^'-3) 

The initial term Hq is given by 

Ho = A-'/' [Hi + J2 Vrniclfrna + fLcOa)]- (5) 
rncr 

In this paper, A is set as 5 and we keep 4500 low-energy 
states for each renormalization step. 

In order to see magnetic properties, we evaluate the 
magnetic susceptibility of /-electron, defined by 



X— 7^ hm 



TrM2j^e-^«/^ Tr^^^e"^"/^' 



, (6) 



where T denotes a logarithmic temperature given by 
r=A^(^~i)/2 ^Yie NRG calculation, H% is the Hamil- 
tonian without H[, and a magnetic moment along the 
2-axis is defined by 
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MB + gsCT/2)fl^^fma + S,,.N- (7) 



Here /ib is the Bohr magneton, Qs is electron g-factor 
given by 3s =2, and 
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(8) 



Note that we consider only the susceptibility for the z- 
component of the magnetic moment, but we have per- 
formed the calculations for all directions of the magnetic 
moment. Then, we have checked that the susceptibilities 
for the magnetic moments along x-, y-, and z-directions 
take the same values due to the cubic symmetry. In the 
following, we show only the results for the susceptibility 
for the magnetic moment along the z-axis. 

In Fig. 2, we show the Curie constant C as a function 
of T. In NRG calculations, it is possible to evaluate pre- 
cisely the Curie constant C by using the relation Tx- We 
note that the finite value of C suggests the emergence of 
of magnetic behavior even for the non-magnetic ground 
state. Here the CEF potential is chosen as W^=— 105 meV 
in order to obtain the small excitation energy. Concern- 
ing the large magnitude of W , we will provide a brief 
comment later. 

First let us consider the case of A=0, in which we find 
C=16/i|/fcB almost irrespective of a temperature, as ob- 
served in Fig. 2. Note that the residual entropy from the 
localized magnetic moment is due to the speciality of the 
present model in which we consider only the hybridiza- 
tion with single Ou band. If we further include iiu and 
bands, the residual entropy will be released due to 
the Kondo effect. 

Now we consider the meaning of "16" in the value of 
C for A=0. In this situation, the ground state is charac- 
terized by S'=3 and L=2t due to the Hund's rules, indi- 
cating 49-fold degeneracy in the ground state. When 
we include the effect of CEF potential, we note that 
the CEF potentials act on the charge, not on the spin. 
Namely, the septet of _L=3 splits into one single and two 
triplets due to the effect of the cubic CEF potentials. 
For < and < a; < 1, the lowest energy state is 
found to be the singlet. Thus, the ground state becomes 
septet characterized by S'=3. As is well known, the Curie 
constant of the susceptibility of the magnetic moment 
M^L + gsS=gjJ is given by C=gj^llJiJ + 1)/(3A:b), 
where gj is the Lande's g-factor. In the present case, 
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J is composed only of electron spins and thus, we set 
J=S=3 and gj=gg=2, leading to C=16^|/(3fcB). Such 
a large value of C at A=:0 is related to the appearance of 
magnetism for A > 0. 

When we increase the value of A, C should be de- 
creased, since the ground state becomes the singlet char- 
acterized by J=0. However, as shown in Fig. 2, the Curie 
constant does not vanish immediately, when we increase 
the value of A. In fact, as observed in the curve for 
A=0.1, C is still constant in the temperature region of 
T > 10"'' due to the smallness of the excitation energy. 
When A is further increased, the temperature dependence 
of the Curie constant becomes significant. For A=0.4, C 
becomes small drastically at a temperature as large as 
T- 10-2. 

In order to compare the numerical results with the ex- 
perimental ones, we plot the inverse susceptibility 1/x in 
Fig. 3. Here we use a large value for CEF potential as 
H^=— 105meV, in order to reproduce the small excitation 
energy. We note that 1/x is shown in the linear scale of 
the temperature, even though the low-temperature pre- 
cise behavior in the logarithmic-scale calculation cannot 
be observed in the linear scale for the temperature. Note 
also that except for A=0, 1/x diverges near T ~ 0, since 
the ground state is singlet and the Curie constant eventu- 
ally goes to zero at low temperatures. In the expression 
of 1/x in the linear scale, even for A=0.3, we find the 
region of Curie- Weiss-like behavior in a temperature re- 
gion of T > 0.01. For A=0.4, it is difficult to observe 
the Curie- Weiss-like behavior in the temperature range 
of T < 0.1. From these results, we conclude that the 
magnetic behavior of Cm02 can be explained even for 
the non-magnetic ground state of Cm^"*" if we appropri- 
ately choose the parameters in the model. 

In the present scenario, the magnetic behavior should 



disappear at low enough temperatures. However, it will 
be difficult to perform low-temperature experiments for 
Cm02. An executable test is to measure the tempera- 
ture dependence of the magnetic moment in the neutron 
diffraction experiment. If the present scenario is correct, 
it is possible to detect the reduction of the moment due 
to the decrease of temperature. 

Finally, we briefly comment on the value of W . In the 
present calculation for x, the value of W may be too large 
in comparison with that for other actinide dioxides.— 
Namely, our result indicates that such a large CEF po- 
tential is required to explain the magnetic behavior in 
Cm02 on the basis of the model calculations. A possi- 
ble detection of the large CEF effect in Cm02 will be 
another test for the present scenario. 

In summary, we have analyzed the magnetic suscepti- 
bility of Cm02 on the basis of the seven-orbital Ander- 
son model by using the NRG technique. We have actu- 
ally found that the magnetic behavior appears even if the 
ground state of Cm ion is non-magnetic, since the excita- 
tion energy between ground and magnetic excited states 
becomes small due to the combined effect of Coulomb in- 
teractions, spin-orbit coupling, and CEF potential. We 
believe that the scenario works to understand the mag- 
netism in Cm02. 
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